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A first-principles approach is presented for the thermoelectricity in molecular junctions formed by 
a single molecule contact. The study investigates the Seebeck coefficient considering the source-drain 
electrodes with distinct temperatures and chemical potentials in a three-terminal geometry junction. 
We compare the Seebeck coefficient in the amino-substituted and unsubstituted butanethiol junction 
and observe interesting thermoelectric properties in the amino-substituted junction. Due to the 
novel states around the Fermi levels introduced by the amino-substitution, the Seebeck coefficient 
could be easily modulated by using gate voltages and biases. When the temperature in one of the 
electrodes is fixed, the Seebeck coefficient varies significantly with the temperature in the other 
electrode, and such dependence could be modulated by varying the gate voltages. As the biases 
increase, richer features in the Seebeck coefficient are observed, which are closely related to the 
transmission functions in the vicinity of the left and right Fermi levels. 



Building electronic circuits from molecules is an in- 
spiring idea [J, 0, 0, Much attention has been de- 
voted to investigating the various transport properties 
that might be applicable in developing new forms of elec- 
tronic and energy-conversion devices, such as electron 
transfer 0, Q , shot noise Qjheat transport [1, 0] , neg- 
ative differential resistance [10[, and gate-controlled ef- 
fects [ll[. Recently, topics on thermo-related transport, 
such as local heating [l^, and thermal transport 

14| , have emerged as new subfields in molecular electron- 
ics. Another important thermo-related property in the 
molecular tunnel junction (m- M-m) is thermoelectricity 

H, [3, [H El, Ei, [111,!!!, The Seebeck coeffi- 
cient, which is related not only to the magnitude but also 
to the slope of the transmission function in the vicinity of 
Fermi levels, can provide more information than current- 
voltage characteristics. The study of thermoelectricity is 
of key importance in the design of novel thermo-related 
electronic and nanoscale energy conversion devices. 

Recent experimental measurements of the Seebeck co- 
efficient for molecular tunnel junctions were conducted 
at zero bias 



21 



l23j . In such cases, the system can be 
described by only a single Fermi level. Similarly, the See- 
beck coefficient for bulk material is also described by a 
single Fermi level. Nevertheless, molecular tunnel junc- 
tions consist of two electrodes as independent electron 
and heat reservoirs. Thus, it is worthwhile extending 
the investigation of the Seebeck coefficients to a system 
with distinct temperatures (Ti^^r)) and chemical poten- 
tials (/iL(i?)) in the left (right) electrode. In this letter, we 
present a theory for two distinct Fermi levels in molecu- 
lar tunnel junction combining a first-principles approach 
for the Seebeck coefficient in the two- and three-terminal 
junctions in nonlinear regime. As an example, we system- 
atically investigate the dependence of the Seebeck coeffi- 
cient on the source-drain biases, gate voltages, and tem- 
peratures in the metal electrodes before and after amino- 
substitution in the butanethiol molecular junction. The 
Seebeck coefficient may provide further insights into the 
physical properties of molecular tunnel junctions. For 



example, whether the Fermi energy is closer to the low- 
est unoccupied molecular orbital (LUMO) or the high- 
est occupied molecular orbital (HOMO) may be locally 
probed via thermoelectricity measurements [l6j . Inter- 
esting features in the Seebeck coefficient are observed 
in the amino-substituted butanethiol junction because 
of the dramatic change in the transmission functions by 
amino-substitution in the vicinity of the left and right 
Fermi levels. 

Alkancthiol [CIl3(CH2)ri-iSH. denoted as C„]-related 
molecules are a good representation of reproducible junc- 
tions that can be fabricated [3, S 01 • It has been estab- 
lished that non-resonant tunneling is the main conduc- 
tion mechanism inasmuch as the Fermi levels of the two 
electrodes lie within the large HOMO-LUMO gap. How- 
ever, functional group substitution may have significant 
effects on the electronic structures of alkanethiols. New 
states around the Fermi levels are produced when -NH2 
is substituted for -H in bridging butanethiol (C4). The 
response of these states to external biases depends on 
the polarity and it leads to asymmetric current-voltage 
characteristics [25[. Because of the dramatic change in 



the transmission function in the vicinity of the Fermi 
levels, the novel characteristics of the Seebeck coefficient 
are observed in the amino-substituted junction. Conse- 
quently, the amino-substitution significantly affects the 
Seebeck coefficient. For example, the Seebeck coefficient 
can change signs by applying gate voltages and biases 
in the amino-substituted junction but not in the un- 
substituted system. The influence of different temper- 
atures between the two electrodes on the Seebeck coeffi- 
cient, controllable by the gate voltages, is significant in 
the amino-substituted junction. The results suggest that 
the thermoelectric molecule devices, such as a molecular 
thermometer, are possible in the future. 

Let us start by considering a single molecule sand- 
wiched between two bulk electrodes applied with a cer- 
tain source-drain bias. The Fermi level in the left /right 
electrodes is determined by filling the conduction band 
with the valence electrons in the bulk Au electrode de- 
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scribed by the jeUium model (r^ « 3). The gate volt- 
age is introduced as a capacitor composed of two par- 
allel circular charged disks separated by a certain dis- 
tance from each other [ll, 27|. The axis of the capacitor 
is perpendicular to the transport direction. One plate 
is placed close to the molecule, while the other plate, 
placed far away from the molecule, is set to be the zero 
reference energy [Inset in FigureHl^a)]. Using the second- 
quantization field-operator technique with the effective 
single-particle wave functions calculated self-consiste ntly 
in density functional theory, the current is given by [2q |: 



\E)-fk{^L,n)r^{E)l{l) 



where the transmission function of electron with energy 
E incident from the left (right) electrode is: 



dR / dK 



\^EE 



(r,K||), (2) 



where Z^^/ 



\I>^^^''(r, K||) is the single-particle wave function (de- 



-V[vE'y**^^„andi, 



J = L, R. 



tailed theory can be found in Ref. i26|,i28i]) incident from 



the left (right) electrode with energy E and component 
of the momentum K|| parallel to the electrode surface, 
and c?R represents an element of the electrode surface. 
The stationary wave function ^'^'•^■'(r, K||) can be cal- 
culated by solving the Lippmann-Schwinger equation it- 
eratively to self-consistency 28, 2^, S^l- The exchange- 
correlation potential is included in density-functional for- 
malism by using the local-density approximation 31 1. 
Once the single-particle wave functions are calculated 
sclf-consistently, the transmission function of electron 
with energy E can be calculated using Eq. ([2]). We as- 
sume that the left (right) electrode serves as the electron 
and thermal reservoir with the electron population de- 
scribed by the Fermi-Dirac distribution function, /^^^' = 
1/ (exp [{E - fJ-L{R)) /kBTL(R)) + l), where iil(r) and 
2l(_r) are the chemical potential and the temperature 
in the left (right) electrode, respectively, and ks is the 
Boltzmann constant. 



This research considers the extra current induced by 
an additional infinitesimal temperature (AT) and volt- 
age (Ay) distributed symmetrically across the molecular 
junction: 



A/ = I{iil.Tl 



AT AT, 

— ■,Hr,Tr — ) + I(pL 
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Tl; fJ-R ~ Tfl) - 2I{pl,Tl; (J-r, Tr), (3) 



The Seebeck coefficient (defined as S* = ■^) is ob- 
tained by letting AI ~ 0. We expand the Fermi-Dirac 
distribution function to the first order in AT and AV 
and obtain 



S 
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where 

i^f = - I dE(E- 



ML(i?)J 



dE 



(4) 



iE), (5) 



and t{E) = t^{E) = t^{E), a direct consequence of the 
time-reversal symmetry. 

The research explores the dependence of the Seebeck 
coefficient on the gate voltages, temperatures in the elec- 
trodes, and the source-drain biases in both the linear and 
nonlinear response regimes by applying Eq. ([4]). In the 
low-temperature regime where the higher order terms in 
the temperature are disregarded, Eq._(|2)_can_be simpli- 
fied using the Sommcrfeld expansion 
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where the Seebeck coefficient is closely related to the 
transmission function in the vicinity of the left and right 
Fermi levels. 

As the first step in our analysis, we study the See- 
beck coefficient in a three-terminal geometry in the lin- 
ear response regime {Vsd = 0.01 V and /i^ w ^r w 
Ep), where both electrodes have the same temperatures 
{Tl = Tr = T). In this case, the Seebeck coefficient 

1 nE-E,)^T{E)dE 



can be simplified as 5* 



imkr{E)dE 



In 



the low temperature regime, the Seebeck coefficient can 
be further simplified using the Sommerfeld expansion as 
S 



■R^k^T dlnrjE) 



\E=Ef 



SB 111- This equation has 



3e dE 

been applied to the study of several atomic and molec- 
ular systems 



\m 21 



The Seebeck coefficient as a func- 
tion of the gate voltage for various temperatures in the 
amino-substituted and unsubstituted butanethiol junc- 
tion is presented in Fig. [ija). The results show that the 
characteristics of the Seebeck coefficient are sensitive to 
the gate voltages in the amino-substituted butanethiol 
junction. The most striking feature is that the molec- 
ular transistor can be converted from n-type to p-type 
by applying the gate voltages. The Seebeck coefficient 
is close to zero at Vg ~ —2.6 V. As the gate voltage 
further decreases, the sign of the Seebeck coefficient be- 
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FIG. 1: (color online) The Seebeck coefficient S in a three- 
terminal geometry with Vsd = O.Oiy: (a) The S versus 
Vg where Tl ~ Tn = T for the amino-substituted [black 
(thick) lines] and unsubstituted [red (thin) lines] butanethiol 
for T = 50 K (solid line), T = 100 K (dashed line), and 
T = 150 K (dotted line). The inset shows the schematic of 
the three-terminal junction. The gate field is applied in a 
direction perpendicular to direction of charge transport, (b) 
The density of states and the transmission function (inset): 
the left pannels for the amino-substituted butanethiol junc- 
tion at Vg = —3.65, —2.60, and 1.72 V; the right pannels for 
the unsubstituted butanethiol junction at Vg = —1.76, —0.58, 
and 2.92 V. (c) The S versus Tl for an amino-substituted bu- 
tanethiol junction for Vg = —3.38, —2.60, and 1.72 V, where 
Tfl = K. 



comes positive (p-type). For the butanethiol molecular 
junction, the characteristic of the carrier remains n-type 
all the time because the sign of the Seebeck coefficient is 
negative. 

To arrive at the physical reason why the gate volt- 
age can efficiently modulate the Seebeck coefficient, the 
DOSs (transmission functions) are plotted as a function 
of energy for the various gate voltages in Fig. [Hb) (In- 
set of Fig. [Ijb)). We observe that the positive (nega- 
tive) gate voltage shifts the LUMO peak towards higher 
(lower) energies. At Vg = —2.6 V, the peak position 
of the LUMO and transmission function align with the 
Fermi levels, implying that d\n.T{E) / dE\E=EF ~ 0- 
Hence, the Seebeck coefficient is close to zero at the 
gate voltage around Vq = —2.6 V. When the gate volt- 
age is tuned at around Vg = 1.72 V, the Seebeck co- 
efficient is negative because d\nT{E)/dE\E=EF > 0. 
Conversely, the Seebeck coefficient is positive because 
d\nT(E)/dE\E=Ei. < at Vg = -3.65 V. Thus, the 
characteristic of the carrier type of a certain molecular 
junction can be converted from n-type (closer to LUMO) 
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FIG. 2: (color online) (a)The Seebeck coefficient S as a func- 
tion of source-drain biases in a two-terminal geometry for 
amino-substituted (the main graph) and unsubstituted (the 
inset in the upper right corner) butanethiol for T = 50 K 
(solid lines), T = 100 K (dashed lines), and T = 150 K (dotted 
lines) . (b) The density of states and the transmission function 
(inset) for various source-drain biases {Vsd = —0.9, —0.5, 
—0.3, 0.2, and 0.8 V) in the amino-substituted butanethiol 
junction. 



to p-type (closer to HOMO) by tuning the gate voltage. 
In the unsubstituted butanethiol junction as shown in 
the right plane of Fig. \V[h) , it is observed that the elec- 
tron transmission function is always small because the 
location of the Fermi levels lies within the large HOMO- 
LUMO gap. We also note that the Seebeck coefficient 
has a negative value due to d\nT{E) /dE\E=EF > 0- 
When the gate voltages are further decreased, the ab- 
solute value of the Seebeck coefficient becomes bigger 
even in the unsubstituted butanethiol junction because 
the negative gate voltage shifts the LUMO peak towards 
the lower energy region. 

The Seebeck coefficient is relevant to the temperatures 
of the electrodes. This property may be applied to the 
design of a molecular thermometer. To show this, we in- 
vestigate the Seebeck coefficient of the amino-substituted 
butanethiol junction at Vsd = 0.01 V as a function 
of temperature of the left electrode (Tl) while keeping 
Tfj = K as shown in Fig. W[,c). The results show that 
the dependence of the Seebeck coefficient on Tl is linear 
at low temperatures. In this regime, the Seebeck coeffi- 
cient can be well described by Eq. where Tr = K. 
As the temperature becomes large, the approxima- 
tion of Eq. ([6|) turns out to be inappropriate, and the 
Seebeck coefficient shows nonlinear behavior. We fur- 
ther observe that the sensitivity of the Seebeck coeffi- 
cient versus Te can be amplified by applying gate volt- 
ages. At Vg — —2.6 V, the Seebeck coefficient has a 
very small value insensitive to Te because the peak of 
the transmission function lies between two Fermi levels. 
When the gate voltage is tuned to Vq = —3.38 V, the 
Seebeck coefficient can be enhanced to around 38 /.tV/K 
at Te = 300 K. The other interesting phenomenon ob- 
served is the possibility to change the sign of the Seebeck 
coefficient by applying the gate voltage. When the gate 
voltage is tuned to Vg = 1-72 V, the Seebeck coefficient 
becomes around -42 ^V/K at Tl = 300 K. The results 
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show that the amino-substituted butanethiol may be an 
effective thermoelectric material applicable to the design 
of molecular thermoelectric devices such as a thermome- 
ter. 

The project further investigates the Seebeck coefficient 
in the nonlinear regime. The Seebeck coefficient of the 
amino-substituted (unsubstituted) butanethiol molecular 
junctions as a function of Vsd is plotted in Fig. [Ifb) (In- 
set of Fig. IHb)). At large Vsd, the difference between 
the left and right chemical potentials becomes significant. 
Thus, the transmission functions in the vicinity of both 
the left and right Fermi levels have important contribu- 
tion to the Seebeck coefficient. In Fig. [IJb) we plot the 
DOSs and transmission functions as functions of the en- 
ergy for various biases. For Vsd > 0, the states be- 
tween the left and right Fermi levels are developed into 
a resonant peak similar to what is found in the elon- 
gated silicon point contact [s^. At large Vsd, the trans- 
mission functions around the left and right Fermi lev- 
els are equally important to the Seebeck coefficient (see 
Eq. (O). When Tl = Tr ^ T, we explain the Seebeck 
coefficient in Fig. by considering the source-drain 
bias Vsd at 0.8, — 0.3, and —0.5 V, respectively. For 
these biases, the contribution to the Seebeck coefficient 
is dominated by the transmission function in the vicin- 
ity of the left Fermi level. Thus, Eq. ([6]) can be fur- 
ther simplified as 5 = _ tt fc^T ^''^'^^) [^^^^ . The See- 
beck coefficient is (negative; zero; positive) at Vsd = 
(0.8; - 0.3; - 0.5) V because dlnT{E) / dE\E=f,^ is 
(>;w; < 0). We also observe that there are more ze- 
roes in the Seebeck coefficient as a function of Vsd in 
Fig. WL.a)- For example at Vsd = 0.2 or -0.9 V, the 
peak position of the transmission function is located in 
the middle of the left and right Fermi levels such that 
dT{E)/dE\E=f,i^ « -dT{E)/dE\E=^n. Consequently, 
the Seebeck coefficient is close to zero at Vsd ~ 0.2 and 
—0.9 V according to Eq. ([6]). In the unsubstituted bu- 
tanethiol junction, the Seebeck coefficient as a function of 
Vsd is shown in the inset of Fig.[2)Ja). The results show 
that the Seebeck coefficient remains a negative value in 
the whole bias regime owing to the fact that the Fermi 
levels are located between the large HOMO-LUMO gap. 

In conclusion, the study investigates the thermoelec- 
tricity in the molecular junction in both linear and non- 
linear regimes. The Seebeck coefficients are studied us- 
ing first-principles calculations. The general properties 
of the Seebeck effects can be very different for the un- 
substituted and amino-substituted butanethiol junction 
in the two-terminal and three-terminal molecular geome- 
tries. The research illustrates that the gate field is able to 
modulate and optimize the Seebeck coefficient. Another 
interesting phenomenon is the possibility to change the 
signs of the Seebeck coefficient by applying the gate volt- 
ages and biases in amino-substituted butanethiol junc- 
tion. It is observed that the Seebeck coefficient is relevant 
to the temperatures of the electrodes that may be applied 
to the design of a molecular thermometer, and its sensi- 



bility can be controlled by gate voltages. We also extend 
the investigation of the Seebeck coefficient to molecular 
tunnel junction at finite biases. As the biases increase, 
richer features in the Seebeck coefficient are observed, 
which are closely related to the transmission functions in 
the vicinity of the left and right Fermi levels. All results 
show that the molecular tunnel junction based on alka- 
nethiols may be a promising candidate for the design of 
novel thermoelectric devices in the future. 

because d\nT{E)/dE |_E=p^ is 
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